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We investigate experimentally and theoretically thin layers of colloid particles held adjacent to 
a solid substrate by gravity. Epifluorescence, confocal, and holographic microscopy, combined with 
Monte Carlo and hydrodynamic simulations, are applied to infer the height distribution function of 
particles above the surface, and their diffusion coefficient parallel to it. As the particle area fraction 
is increased, the height distribution becomes bimodal, indicating the formation of a distinct second 
layer. In our theory we treat the suspension as a series of weakly coupled quasi-two-dimensional 
layers in equilibrium with respect to particle exchange. We experimentally, numerically, and theo¬ 
retically study the changing occupancies of the layers as the area fraction is increased. The decrease 
of the particle diffusion coefficient with concentration is found to be weakened by the layering. 

We demonstrate that particle polydispersity strongly affects the properties of the sedimented layer, 
because of particle size segregation due to gravity. 


I. INTRODUCTION 

Being relevant to a wide range of practical scenarios, 
the behavior of colloid suspensions near solid surfaces has 
been thoroughly studied over the years. This research ef¬ 
fort consists of several bodies of work, for each of which 
we can give only a few representative references. The 
first category of papers concerns the disruption of the 
structural isotropy of a three-dimensional (3D) fluid sus¬ 
pension by the surface, e.g., the formation of a layered 
structure decaying away from the surface under equilib¬ 
rium [Il,i and nonequilibrium Q conditions. Another 
category addresses the effect of the anisotropic geometry 
on particle dynamics near a single planar surface — for 
isolated particles [MU, particle pairs iJMl , and a 
3D suspension adjacent to a surface 

Regarding quasi-two-dimensional (quasi-2D) layers of 
particles, most studies have considered the confinement 
of suspensions between two rigid surfaces. This research 
addressed structural properties of such confined suspen¬ 
sions [20l - [^ , and the dynamics of single particles |25l- 
[^ . particle pairs [^ . and concentrated quasi-2D 
suspensions [3(M^ . Another type of quasi-2D suspen¬ 
sions has also been studied, where a particle layer is con¬ 
fined to a fluid interface [ssl - fs^ . 

In cases where the surface attracts the particles and the 
suspension is sufficiently dilute, the system can contain a 
single layer of surface-associated particles in contact with 
a practically particle-free solvent 0. A single layer can 
also form as a result of gravitational settling of particles 
toward a horizontal wall. This scenario is studied in the 
present work. 

Sedimented colloidal particles undergo random Brow¬ 
nian displacements, which results in diffusive broadening 
of the fluctuating particle layer. The width of the parti¬ 
cle height distribution above the bottom surface is char¬ 
acterized by the sedimentation length I, i.e., the height 


at which the gravitational energy of a particle equals its 
thermal energy. The dynamics and height distribution 
of individual sedimented particles above the bottom sur¬ 
face were studied in Refs. using total internal re¬ 
flection microscopy. Particle monolayers at higher den¬ 
sities were investigated experimentally for a system in 
which the sedimentation length is much smaller than the 
particle diameter [^. It was shown that at high area 
fractions the suspension can assemble into quasi-2D col¬ 
loidal crystals, but formation of a nonuniform vertical 
microstructure was not observed, because of the small 
sedimentation length. 

Here we are interested in the structure and dynamics 
of a surface-associated layer for which the sedimentation 
length is comparable to the particle diameter. We focus 
on the effects of the suspension concentration on the sta¬ 
tistical height distribution of particles and their diffusion 
coefficient. Unlike the quasi-2D suspensions confined be¬ 
tween two surfaces or adsorbed at a fluid interface (which 
restricts particle configurations and motions in two direc¬ 
tions) in the present system no constraints are imposed 
on the distance between the particles and the single wall. 
Thus, at sufficiently high area fractions, particles form a 
nontrivial stratified microstructure. This microstructure 
and its effect on particle dynamics are analyzed in our 
paper. 

The article is organized as follows. Section HIl describes 
the experimental methods used to prepare the system, 
image the particles, and analyze the extracted data. In 
Sec.lm] we describe the theoretical background and nu¬ 
merical methods used to perform the simulations. In 
Sec. II VI we present the results concerning the equilibrium 
structure of the quasi-2D suspension observed in planes 
parallel to the bottom surface (the quasi-2D radial dis¬ 
tribution) and in the direction perpendicular to it (the 
height distribution). Section IVl addresses the diffusion of 
particles parallel to the surface, as affected by the surface 
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Figure 1. (a) Images of fluorescent 1.5 /rm-diameter silica 

spheres suspended in water, taken after the particles sed¬ 
imented to create a quasi-2D suspension at area fraction 
(/> = 0.49. Large (small) image corresponds to a typical image 
of the hrst (second) layer. Scale bar = 5 pm. (b) Schematic 
view of the system and its parameters. 


proximity. We discuss our findings in Sec. IVII 


II. EXPERIMENTAL METHODS 

A. Quasi-2D system of sedimented Brownian 
spheres 

Quasi-2D colloidal layers are created by placing a 
suspension of colloidal silica spheres in a glass sample 
cell ~ 150 pm high. The particles are then allowed to 
sediment and equilibrate for 30 minutes at a tempera¬ 
ture of approximately 24°C before measurements start 
(Fig. [T]). We use green fluorescent monodisperse, nega¬ 
tively charged silica particles (Kisker Biotech, PSTG1.5 
Lot #GK0090642T) with diameter d = 1.50 ± 0.15 pm, 
and mass density go = 2.0 g/cvo?. Monolayers of area 
fraction 0 < (j) < 0.62 are prepared by diluting the 
original suspension with double distilled water (DDW, 
18 Mfl), without and with the addition of salt at a con¬ 
centration [KGl] = 0.01 M. The sample walls are cleaned 
and slightly charged by plasma etching to avoid particle 
attachment to the bottom wall of the cell. We observe 
that the aqueous medium above the colloidal monolayer 
is free of colloids. Since the particles are floating right 
above the bottom wall, we can treat the upper wall as a 
distant boundary. 



Figure 2. Holographic imaging. Height probability distri¬ 
bution of a single sphere (in salt-free water) shifted to the 
maximum value and fitted to the Boltzmann probability dis¬ 
tribution with the particle-wall potential ©• 

B. Imaging techniques 

Particle position and motion in the x-y plane, perpen¬ 
dicular to the optical axis, are observed using epifluo- 
rescence microscopy (Olympus 1X71). Images are cap¬ 
tured at a rate of 70 fps by a CMOS camera (Gazelle, 
Point Grey Research). We use in-line holographic mi¬ 
croscopy to image the dynamics of particles in three di¬ 
mensions in dilute samples [s^. This imaging technique 
uses a collimated coherent light source (DPSS, Coherent, 
A = 532nm) to illuminate a sample mounted on a micro¬ 
scope. The light scattered from the sample interferes 
with the light passing through it, to form a hologram 
in the image plane. We reconstruct the light field pass¬ 
ing through the sample by Rayleigh-Sommerfeld back- 
propagation and extract from it the particle location in 
three dimensions [38l - l40l| . For holographic imaging mea¬ 
surements we use non-fluorescent silica particles with the 
same diameter (d = 1.50 ± 0.08 /rm. Polysciences Inc.). 
Additional details of the setup and measurement meth¬ 
ods can be found elsewhere |38l| . 

We use confocal imaging to monitor particle posi¬ 
tions in a dense layer in three dimensions. Our spin¬ 
ning disc confocal imaging system (Andor, Revolution 
XD) includes a Yokogawa (CSU-Xl) spinning disc, and 
an Andor (iXon 897) EM-CCD camera. An objective 
lens (Olympus, x60, NA=1.1, water immersion) mounted 
on a piezoelectric scanner (Physik Instrumente, Pifoc P- 
721.LLQ) is used to scan the sample in the z axis, with 
a step size of 100 nm. 

C. Height calibration 

A suspended tracer particle is subject to electrostatic 
and gravitational forces in addition to thermal fluctua¬ 
tions, affecting its height distribution The particle 
potential energy can be described as 

U = mgz + (1) 
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where z is the vertical position of the tracer, g is the 
gravitational acceleration, 

m ( 2 ) 

is the buoyant mass of the tracer (Apo is the mass den¬ 
sity difference between silica and water), A is the Debye 
screening length, and the amplitude B depends on A and 
the surface charges of both particle and glass surfaces. 
The corresponding probability distribution of the parti¬ 
cle height z is 

p{z) = ^3) 


where k^T is the thermal energy and Z~^ is the normal¬ 
ization constant. 

The height distribution of a single particle above the 
sample’s bottom was obtained from very dilute sus¬ 
pensions, using in-line holographic imaging [38l - l40l| (see 
Fig. 12). Our holographic measurements provide values of 
relative particle positions, but not the absolute particle 
heights with respect to the bottom wall. We thus set the 
peak position to z = 0 and focus on the height relative 
to this reference plane. The exponential decay on the 
right side of the probability-density peak is governed by 
a decay length. 


ksT 

mg 


( 4 ) 


(the sedimentation length), resulting from the com¬ 
petition between gravity and thermal forces. The 
exponential-decay length determined from the holo¬ 
graphic measurements agrees well with the calculated 
sedimentation length Q , without any fitting parameters 
(see Fig. [2]). The electrostatic term of the probability- 
density, which controls the steep rise of the probability, 
affects mostly the peak position rather than its shape. 
Since we shifted the peak position to 2 ; = 0, the fitting 
of the entire probability-density using Eqs. m and © 
was insensitive to the value of B. Reasonable fits were 
obtained for A in the range of A ^ 40 — 70 nm. Better 
estimations of B and A are given in Sec. IIV A 1[ using 
mobility measurements. 

The applicability of the holographic imaging is limited 
to low-density suspensions, whereas the confocal imaging 
can be also used at higher concentrations. On the other 
hand, confocal height measurements suffer from spherical 
aberrations due to multiple changes in refractive index 
in the imaging path. This leads to a systematic error 
in measuring z, which can be eliminated by proper cal¬ 
ibration. We calibrate the confocal measurement of the 
relative vertical particle positions by requiring the expo¬ 
nential decay of the height distribution to agree with the 
known, and verified, value of 1. 

In Fig. [3Ka) we show the particle-height distribution 
p{z) at (j) < 0.003 for two different particle sizes (do = 
1.0,1.5 pm). The distributions are shifted so that the 
highest probability is located at z = 0. Scaling the loga¬ 
rithm of the distributions by dj) [inset of Fig. EKa)] shows 
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Figure 3. (a) Confocal imaging. Height probability distri¬ 
bution of silica particles at (j) < 0.003 with diameters of 
do = 1.5 pm (blue) and do = 1.0 pm (black). Inset: Log¬ 
arithm of the probability distributions scaled by dp in units 
of 10^ pm~^-, as expected, the two curves have approximately 
the same slope, which is used to calibrate the confocal height 
measurements, (b) Height probability distribution of silica 
particles with diameter do = 1.5 pm extracted from holo¬ 
graphic imaging (green circles, see Fig.[2ll and confocal imag¬ 
ing (blue line). The suspensions in both figures were with no 
added salt, [KCl] = OM. 


that the normalized decay constants for the two parti¬ 
cle sizes have approximately the same value, from which 
we calibrate the confocal microscope’s height measure¬ 
ments. In Fig. [1Kb), the height distributions extracted 
by the two methods (holographic and confocal imaging) 
are overlaid. This figure emphasizes the higher accuracy 
of holographic imaging over confocal imaging, especially 
around 2 : = 0, where the increase in distribution should 
be very steep [1, [s^ . The difference between the curves 
can also be attributed to polydispersity, since the holo¬ 
graphic imaging is a single-particle measurement while 
the confocal imaging is a multiple-particle measurement, 
and its corresponding curve represents an average over 
~40 particles. 
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III. NUMERICAL METHODS 
A. The system 

1. Particles and their interactions 

Silica particles are modeled as Brownian hard spheres 
with or without electrostatic repulsion (depending on the 
salt concentration), immersed in a fluid of viscosity rj. 
The bottom wall is treated as an infinite hard planar 
surface. Creeping-flow conditions and no slip boundary 
conditions at the particle surfaces and at the wall are 
assumed. 

In a salt solution with [KCl] = 0.01 M, the Debye 
length is only about 5 nm, and therefore electrostatic 
interactions are screened out. The particles thus interact 
only via infinite hard-core particle-particle and particle- 
wall potentials and the gravity potential mgz, and no 
other potential forces are involved. The strength of the 
gravity force is described by the sedimentation length ((T|) . 

In addition to the hard-core repulsion, in DDW with 
no added salt ([KCl] = 0 M) particles are assumed to also 
interact via particle-wall and particle-particle Debye- 
Hiickel potentials. 


V{z) = 

(5) 

1 

T 

II 

(6) 


where A is the Debye screening length, B and B' are the 
potential amplitudes, and r is the distance between the 
particle centers. The consideration of Debye-Hiickel po¬ 
tentials in the salt-free case is based on our experimental 
measurement A ^ 60 nm. A finite Debye screening length 
in DDW stems from the presence of residual ions in the 
solution (4ll| . 


where 

mo='^dll\QQ. (9) 

The area fraction (j) based on the average particle diam¬ 
eter do is 

(j)=\TTndl, ( 10 ) 

where n is the number of particles per unit area. Since 
the particles are free to move in the z direction, the area 
fraction (jj can exceed 1. 

3. System parameters 

The simulations were carried out for the following sys¬ 
tem parameters: For the dimensionless sedimentation 
length ([5|) we use the value 

^=0.158, (11) 

calculated from the particle size and density. Based on 
the comparison between the calculated and measured val¬ 
ues of the equilibrium average of the lateral self-diffusion 
coefficient for isolated particles in DDW, we estimate 
that the Debye length and the amplitude of particle-wall 
electrostatic repulsion are 

A/d = 0.03, 7:^ = 10. (12) 

ks-i 

These values are used for salt-free suspensions at all sus¬ 
pension concentrations. Assuming that the charge densi¬ 
ties of the particle and wall surfaces are similar, we take 

B' = B/2, (13) 


2. Suspension polydispersity 


To determine the effects of the suspension polydisper¬ 
sity on the near-wall microstructure and dynamics, we 
have performed numerical simulations for a hard-sphere 
(HS) system with a Gaussian distribution of particle di¬ 
ameters, 


p{d) 


1 

(27ra2)i/2 


exp 


(d- do)^' 


(7) 


where d and do are the actual and average particle diam¬ 
eters, and (7 is the standard deviation. All the particles 
have the same mass density po! hence, particles of dif¬ 
ferent sizes have different buoyant masses and different 
sedimentation lengths (|3]). The dimensionless sedimen¬ 
tation length based on the average particle diameter do, 
is defined as 

^ ksT 
do mogdo ’ 


for the interparticle-potential amplitude, as follows from 
the Derjaguin approximation [d^. 

The simulations were performed in the range of area 
fractions </> < 1.2. For polydisperse HS systems the cal¬ 
culations were carried out for cr/do = 0.10, 0.15, 0.20, 
and 0.25 (we estimate that 0.10 < cr/do < 0.15 for the 
silica particles used in the experiments). For particles 
interacting via the Debye-Hiickel potentials ([5]) and ([6]) 
only monodisperse suspensions were considered. 


B. Evaluation of the equilibrium distribution 

1. Low density limit 

For monodisperse suspensions at low particle concen¬ 
trations, the equilibrium particle distribution p{z) is 
given by the normalized Boltzmann factor (|31). To de¬ 
termine the particle distribution for a dilute polydisperse 
suspension, the particle-size-dependent Boltzmann factor 
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for individual particles, pi{z;d), is convoluted with the 
particle-size distribution 

piz)=f ddp{d)pi{z-d), (14) 

^0 

For a HS system 

pi{z] d) = - d/2), (15) 

according to equations 0-0, where 9{x) is the Heav¬ 
iside step function, and the sedimentation length I is 
particle-size dependent due to the variation of particle 
mass. 


2. Monte-Carlo simulations 

To determine the equilibrium microstructure of a sed¬ 
imented suspension at finite particle area fractions, equi¬ 
librium Monte-Carlo (MC) simulations were performed 
for 2D-periodic arrays of spherical particles in 3D space 
(with periodicity in the horizontal directions x and y and 
the box size L). The particles interact via infinite hard¬ 
core repulsion and the pair-additive potential 


standard deviation adaptively adjusted to the cur¬ 
rent mean gap between particles. This displace¬ 
ment results in the change of the configuration from 
X to X'. 

(b) According to the Metropolis detailed balance con¬ 
dition, the new configuration is accepted with the 
probability 

min (1, exp {- [C/(X') - C/(X)] /keT }), (17) 

provided that there is no particle-particle or 
particle-wall overlap. 

To let the system reach an equilibrium state Xi, the 
MC step (a) and (b) is repeated 10^times. The next 
independent equilibrium configuration Xn+i is obtained 
from the previous configuration by performing lO^iV 
MC steps. The particle height distribution p(z) and other 
equilibrium quantities are obtained by averaging over 200 
independent configurations Xi. 

C. Hydrodynamics and self-diffusion 

1. Low density limit 


N N N N 

U (X) = ruipzi + ^ ( 16 ) 

i=l i=l j/i 


In the absence of a wall, the self-diffusion coefficient of 
an isolated solid sphere with diameter do is given by the 
Stokes-Einstein expression 


which includes the gravity term and particle-wall and 
particle-particle screened electrostatic potentials m and 
®. Here X = (ri,..., r^r) is the particle configuration 
(with denoting the position of particle i), Zi is the 
vertical coordinate of particle i, and = jr^ — rjj is the 
relative particle distance. 

A purely HS system with V = V = Q was modeled 
for monodisperse particles and for polydisperse particles 
with the Gaussian size distribution 0. For systems with 
nonzero electrostatic repulsion only monodisperse parti¬ 
cles were considered. 

The initial configuration was prepared by placing N = 
400 particles randomly in a vertical cuboid box with the 
square base L and the height lOL. The size L of the 
2D-periodic cell was determined to obtain the required 
area fraction </> of the sedimented particle layer. The sus¬ 
pension was allowed to sediment by following the MC 
random-walk dynamics in the configurational space X 
( 4 ^ (as described below). After the equilibrium state 
was reached, suspension properties were obtained by av¬ 
eraging the quantities of interest over at least 200 inde¬ 
pendent configurations. 

Our adaptive simulation procedure was performed by 
repeating the MC steps defined as follows: 

(a) A randomly selected particle i is given a small ran¬ 
dom displacement, —>■ r' = Ti + A, where A is 

chosen from a 3D Gaussian distribution with the 


Smjdo 


(18) 


The self-diffusion coefficient D(z) of a sphere with di¬ 
ameter d at a distance 2 from the wall is smaller by a 
factor 


= (19) 

where the normalized mobility coefficient ^|| depends on 
the dimensionless particle position z/d and no other pa¬ 
rameters. Relation dni) refers to the lateral component of 
the self-diffusion coefficient (parallel to the wall), which 
was measured in our experiments. However, an analo¬ 
gous expression also holds for the normal component. 

For monodisperse particles in the dilute-suspension 
limit, the effective self-diffusion coefficient Ds averaged 
across the suspension layer is obtained by integrating m 
with the Boltzmann distribution 

TT = / <^zp{z)lJ-\\{z/d). (20) 

^0 Jd /2 

For a polydisperse suspension, an additional average over 
the particle-size distribution o is needed, 

n r°° r°° 

-T^= ddp{d) dzpi{z;d)p.i\{z/d). (21) 
^0 Jo Jd/2 
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The mobility coefficient yL\\{z/d) was evaluated with 
high accuracy using the Hydromultipole algorithm for 
a particle near a single wall P. The integrals in Eqs. (12(111 
and (ED) were performed numerically using the Gauss 
method, with p,||(z/d) calculated by a series expansion. 


2. Computations for larger densities 

The effective self-diffusion coefficient for suspensions 
at higher concentrations was evaluated using a periodic 
version of the Cartesian-representation al gori thm for 
a system of particles in a parallel-wall channel [^,1^. In 
our approach, periodic boundary conditions in the lateral 
directions are incorporated by splitting the flow reflected 
by the particles into a short-range near-field contribution 
and a long-range asymptotic Hele-Shaw component. The 
near-field contribution is summed explicitly over neigh¬ 
boring periodic cells, and the Hele-Shaw component is 
evaluated using Ewald summation method for a 2D har¬ 
monic potential [H, ■ 

The one-wall results were derived from the two-wall 
calculations using an asymptotic procedure based on the 
observation that in the particle-free part of the channel 
the velocity field tends to a combination of a plug flow 
and a shear flow. All other flow components decay expo¬ 
nentially with the distance from the particle layer. The 
one-wall results are obtained by eliminating the shear 
flow and retaining only the plug flow generated by hy¬ 
drodynamic forces induced on the particles [d^. The 
calculations were performed for the distance to the up¬ 
per virtual wall H = lOdoj which is sufficient to obtain 
highly accurate one-wall results. 

The self-diffusion coefficient is determined by averag¬ 
ing the trace of the lateral translational-translational 
A^-particle mobility, evaluated using Hydromultipole 
codes based on the above alg orithm, with the multipole 
truncation order L = 2 |49l |. The averaging was per¬ 
formed over equilibrium configurations of iV = 400 par¬ 
ticles in a 2D-periodically replicated simulation cell. In¬ 
dependent equilibrium configurations were constructed 
using the MC technique described in Sec. IHIB2I 


IV. STRUCTURE OF THE QUASI-2D 
SUSPENSIONS 


A. Experimental results 


A typical image of our quasi-2D colloidal suspension is 
shown in Fig. [TJ For each area fraction </> and salt con¬ 
centration, the suspension can be characterized by the 
structure in the x-y plane (parallel to the cell floor) and 
the density profile in the 2 direction (perpendicular to 
the floor). In this section we discuss results of our mea¬ 
surements of the microstructure of a sedimented particle 
layer. 


1. Mean particle height at low area fractions 


As mentioned in Sec. Ill Cl our imaging techniques 
do not yield absolute particle heights. To estimate the 
mean particle distance z from the bottom wall (the mean 
height) in a dilute suspension layer, we observe particle 
dynamics in the horizontal directions, and compare mea¬ 
surement results with theoretical calculations of the effect 
of the wall on the lateral particle diffusion. Using fluo¬ 
rescence imaging, we determine the projection of particle 
trajectories onto the x-y plane, r|| (t), and extract the ef¬ 
fective self-diffusion coefficient. 


D, = (Ar2(r))/(4r), (22) 


where r is the time interval. The position-dependent 
diffusivity D(z) in the x-y plane of a single particle near 
a planar wall is given W the following expansion in the 
particle-wall distance [J, Q , 


Djz) ^ 9d 1 fdy 
Do 32z 64 \z) 

45 /dV 1 /dy 
~ 40% ’ 


(23) 


where z = 0 is the wall position. The expansion (051) is 
accurate within 5% to 1% as z increases from 0.5Id up to 
d/2 + 21, in the range where sedimented particles spend 
most of the time in a low-density suspension under equi¬ 
librium conditions. Here I ~ O.Ibd is the sedimentation 
length (|3|). 

From expression (|23)) and Dg extracted according to 
(051) . we can calculate the suspension’s mean distance 
from the wall (where z in (1231) is replaced by a mean value 
(z)). This calculation holds in the limit ^ 0, where 

there are no particle-particle interactions. We measure 
Dg from the particle trajectories, r||(t), in extremely low 
area fraction solution, (f) < 0.003 (in salt-free water), and 
obtain a mean distance from the wall (z) = I.l ±0.1 /rm, 
corresponding to a mean gap e = z — d/2 of 0.3-0.4/rm 
between the particle surface and the wall. We also ex¬ 
tract (z) for different salt concentrations by extrapolat¬ 
ing Dg (measured at various area fractions) to </> = 0 (see 
Sec. IV Al) . obtaining (z) = 0.95 ± 0.05/rm for [KCl] = 
0.01 M and (z) = 1.11 ± 0.05/rm for [KCl] = OM. The 
latter matches the average height extracted from the dif¬ 
fusion of tracers in the extremely low density suspension. 

For A = 5 nm (added salt), the mean height calculated 
from the Boltzmann distribution ([3]) is dominated by the 
exponential decay due to gravity and is practically in¬ 
dependent of B in the particle-wall potential ([D- For 
B = 0, using the particle mass as determined from Eq. 
([5]) with no fitting parameters, we get (z) = 0.99 /rm, 
in agreement with the diffusivity-based measurements of 
(z). This result confirms that in the added-salt case we 
can neglect the electrostatic repulsion from the wall. For 
the salt-free case, taking A = 50 nm, we obtain (z) = 1.11 
/im for B in the range 5-15 k^T. These values are 
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Figure 4. Radial distribution function g(r) in the x-y plane 
for experiments in (a) salt-free and (b) salt-added ([KCl] = 
0.01 M) water. The distribution g{r) was calculated sepa¬ 
rately in the first and second layers (see Sec. IIV A 41 for the 
layer definition) and combined with appropriate weights. 

consistent with those obtained by fitting the measured 
height distribution to the theoretical expression o (see 
Sec.lITC]). 

Since Eq. (1231) does not include lubrication correction 
for small particle-wall gaps e, it overpredicts D{z) for 
z < O.Sld; however, the accuracy of the approximation 
is sufficient for the purpose of the present estimates. In 
our calculations discussed in Secs. Im] and IVH highly 
accurate Hydromultipole results were used instead of 
the far-held approximation (051) . 


2. Radial distribution in the horizontal plane 

To verify that no crystalline or hexatic structures are 
formed at higher values of the area fraction, we evaluate 
from the experiment the radial distribution function g{r) 
and the full 2D pair distribution g{r, 9) in the x-y plane, 
for both the salt-free and salt-added suspensions. No 
dependence on 9 was found. The radial distribution g{r) 
for several values of the area fraction (j) is shown in Fig. 
SKa) for the salt-free system and inSKb) for the salt-added 
system. 

For inonodisperse hard spheres the hrst peak of g{r) 
should correspond to the diameter of the sphere. Our 
measurements show that the hrst peak is at r = 1.68 ^m 


for suspensions without salt and at r = 1.60 /iin for sus¬ 
pensions with [KCl] = 0.01 M. The difference between 
these two numbers implies that the effective shell around 
the particles in the salt-free samples is around 40-50 nm, 
which provides an estimation for the screening length in 
DDW without the addition of salt. This estimate of A is 
consistent with the other two mentioned above. 


3. Vertical density profile 

The height distributions p{z) of the silica particles at 
diherent area fractions of the sedimented particle layer 
were acquired using confocal imaging and conventional 
image analysis [s^. These distributions for salt-added 
suspensions with [KCl] = 0.01 M are plotted in Fig. [SKa) 
for several values of the area-fraction fi. Since we cannot 
precisely measure the position of the wall, the distribu¬ 
tions are shifted so that their first peak (close to the wall) 
is located at z = 0. These distributions indicate the for¬ 
mation of a second layer of particles for area fractions 
(p > 0.26. The observed center of the second layer is lo¬ 
cated Az ss 0.75/rm above the center of the first layer. 
The layer separation is thus significantly smaller than the 
expected separation Az « d = 1.5 /rm (which is similar 
to the peak separation for the radial distribution). See 
further discussion in Secs. ED and Appendix El 

To highlight the onset of the formation of the second 
layer, we look at the subtraction of the height probability 
distribution of the lowest area fraction from the distribu¬ 
tion of all area fractions, Ap = p — p<^=o.o 54 [Fig- Ela)]- 
Two phenomena are expected when a second layer is 
formed: (i) negative values at z = 0 /rm, corresponding to 
a reduction in the fraction of particles populating the first 
layer, (ii) positive and increasing values at z = 0.75 /iin, 
corresponding to the formation and increasing popula¬ 
tion of the second layer. The values of Ap at z = 0 and 
0.75 /im are plotted in Fig.|6Kb). The two expected phe¬ 
nomena are observed at approximately (p ~ 0.3, indicat¬ 
ing the area fraction above which a second layer becomes 
occupied. At area fractions smaller than 0.3 we still ob¬ 
tain negative values of Ap at z = 0 /rm, and positive 
values at z = 0.75 /rm, however these values are rela¬ 
tively low and can correspond to the broadening of the 
exponential distribution due to increase in cp. 


4- Particle-layer occupation fractions 

For the area fractions at which a clear second peak in 
the particle distribution p is seen in Fig. [SKa) (i.e., for 
(p = 0.48 and 0.54), we ht the area around each peak 
to a Gaussian function and define the point of intersec¬ 
tion between the two Gaussians as the effective boundary 
between the two layers. Figure EKb) shows the two dis¬ 
tributions with the Gaussian fits and our definition of 
that boundary, which turns out to be at a distance of 

















Figure 5. (a) Height probability distribution of the silica col¬ 
loids (in [KCl] = 0.01 M) for increasing area fraction reveals 
the formation of a second layer. Colors correspond to different 
area fractions (as labeled), (b) For the most dense suspen¬ 
sions [<j) = 0.54 and (fi = 0.48 (inset)], the height distribution 
(black solid line) around the two peaks can be fitted to two 
Gaussian functions (blue lines). The intersection of the two 
Gaussians defines an effective boundary (red broken line) be¬ 
tween the first and second layers; occupation percentages are 
indicated. 


0.39 ± 0.04 fim above the peak of the first layer in both 
densities. 

Using this boundary, we evaluated the occupation frac¬ 
tions fi = (tii/cf) of the bottom [i = 1) and top layer 
(i = 2), where (j)i is the area fraction of particles in 
layer i. The results are shown in Fig. EKa) for a sus¬ 
pension in [KCl] = 0.01 M solution as a function of the 
total area fraction (j). As expected, the fraction of parti¬ 
cles populating the second layer grows as the total area 
fraction of the suspension is increased. 

An additional independent measurement of the layer 
occupation fractions is done using epifluorescence mi¬ 
croscopy, which enables us to image the different layers 
separately [see Fig.[T][a)]. The occupation fraction of each 
layer is determined by counting the number of particles 
observed therein. The occupation fractions measured us¬ 
ing the epifluorescence imaging technique are plotted in 
Fig. EKa) along with the results obtained from the confo- 
cal microscopy. The two methods yield similar results. 

Alternatively, we can represent the layer-occupation 
results in terms of the area fractions (j)i and (j32 of the 



Figure 6. (a) The difference between the height probability 
distribution at increasing area fractions and the distribution 
at the lowest area fraction (j) = 0.054, Ap = p — pp=o .054 (in 
[KCl] = 0.01 M). Colors are as in Fig[5l[a). Gray (black) 
dashed line corresponds to z = 0.75 pm {z — 0 pm), (b) 
Values of Ap at z = 0 pm (red squares) and 2 = 0.75 pm (blue 
circles) for all area fractions. Both plots exhibit a change in 
trend at area fraction f ~ 0.3. 


first and second layers [see Fig. Elb)]- Both (f>i and (j )2 
increase as is increased, and (pi seems to saturate at 
p > 0.45. 


B. Numerical simulations 

Here we present results of MC simulations of the equi¬ 
librium microstructure of a HS suspension in the near 
wall region. The HS potential corresponds to the system 
with [KCl] = 0.01 M, for which the electrostatic repulsion 
is negligible. Since the suspension used in the exper¬ 
iments is polydisperse, we consider both monodisperse 
and polydisperse systems. 


1. Near-wall particle distribution 

Figure [H] shows MC results for the suspension density 
profile p{z) for a monodisperse suspension and polydis¬ 
perse suspensions at several area fractions. Similarly to 
the experimental results, the simulations show that there 
is a single layer of sedimented particles at low area frac¬ 
tions p, and a two-layer microstructure at higher area 













9 



4 * 



4 * 

Figure 7. Occupation fractions fi and area fractions cpi of 
suspension layers as a function of the total area fraction 0; 
circles (squares) correspond to the first (second) layer, (a) 
Occupation fractions for experiments in [KCl] = 0.01 M ex¬ 
tracted from the confocal height distribution curves (yellow) 
and from the 2D images (green), showing good agreement be¬ 
tween the two methods, (b) The same data replotted for the 
area fractions cpi and 4>2 of the first and second layers. 


fractions. (Development of a third layer for (j) > 0.9 
is also noticeable in the region z/dg > 2.) Suspension 
polydispersity results in broadening of the peaks of the 
particle distribution. 

A direct comparison between the experimental and 
simulation results is presented in Fig. [5] for two values of 
the area fraction (j). At low area fractions [Fig. IHKa)] the 
agreement between the experiments and simulations is 
good. (The standard deviation of the particle-size distri¬ 
bution for which the simulations match the experimental 
data, a/do ~ 0.25, is larger than the estimated stan¬ 
dard deviation 0.1 < a/do < 0.15 based on the manu¬ 
facturer’s specifications; the additional spread of the ex¬ 
perimentally observed peak can be attributed to random 
errors of the particle height evaluation from the confocal- 
microscopy images.) 

A comparison of the numerical and experimental re¬ 
sults at a higher area fraction, as shown in Fig. HKb) 
[also see Figs. |S] and [H], reveals that (j) the experimen¬ 
tally observed second maximum of the density distribu¬ 
tion develops at lower area fractions than the correspond¬ 
ing maximum in the numerical simulations; {ii) the ex¬ 
perimental second peak is narrower, and its position is 



z/do 



z/do 


Figure 8. Particle-wall distribution function for (a) monodis- 
perse suspension; (b) polydisperse suspension with standard 
deviation of particle diameter a/do = 0.15. Simulation results 
for area fraction </> = 0 (solid line), 0.3 (dashed), 0.6 (dotted), 
0.9 (dot-dashed), 1.2 (long-dashed). The insets show the de¬ 
viation Ap — p — p 4 ,^o from the low-density distribution. 


shifted towards the wall. In contrast, the plots of the 
excess distribution Ap with respect to the low-density 
limit, shown in Fig. UKa) and the insets of Fig. [51 indi¬ 
cate that the onset of the formation of the second layer 
occurs at approximately the same area fraction accord¬ 
ing to the simulations and experiments. Moreover, the 
measured and calculated occupation fractions of the lay¬ 
ers are similar for all area fractions, as depicted in Fig. 
fTOl' a). A possible source of the observed discrepancies 
between the experimental and numerical results for the 
particle distribution p{z) is described in Appendix [^ It 
also provides a plausible explanation of the fact that the 
agreement between the experiments and MC simulations 
for the layer occupation fractions fi is quite good in spite 
of the discrepancies for p{z). 


2. Polydispersity effects 

The results in Fig. fTUI a) show that the occupation frac¬ 
tion of the first two layers is relatively insensitive to the 
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Figure 9. Particle-wall distribution function for area frac¬ 
tions (a) (j) = 0.054 and (b) 0.54. Experimental results (solid 
circles); simulation results for standard deviation of particle 
diameter cr/do = 0 (solid line), 0.1 (dashed), 0.15 (dotted), 
0.20 (dashed-dotted), and 0.25 (long-dashed). 


suspension polydispersity; in contrast, the occupation 
fraction of the third layer strongly increases with the 
standard deviation of particle diameter. This increase 
stems from the presence of smaller (lighter) particles in 
polydisperse systems: smaller particles tend to migrate 
into the top layer, as evident from Fig. irUT b). For di¬ 
lute suspensions, the particle-size segregation results in 
variation of the slope of log p[z) with the distance from 
the wall, as illustrated in Fig. [TTJ We estimate that this 
variation causes an approximately 20% uncertainty of 
the calibration of the confocal height measurements de¬ 
scribed in Sec. Ecl 


C. A quasi-2D model of the equilibrium layered 
microstructure 

Here we present a semi-quantitative theoretical model 
for evaluating the occupation fractions fi of the particle 
layers in a sedimented colloidal suspension. Our theory is 
based on the assumption that the suspension inicrostruc- 
ture can be approximated as a collection of weakly cou- 


Figure 10. (a) Occupation fraction fi and (b) normalized av¬ 
erage particle diameter di in the first particle layer (black), 
second layer (blue), and third layer (green), vs the total area 
fraction 0. The results for a monodisperse system (solid 
lines) and polydisperse systems with standard deviation of 
particle diameter cr/do = 0.1 (dashed), 0.15 (dotted) 0.20 
(dash-dotted), and 0.25 (long-dashed). The symbols repre¬ 
sent experimental results from confocal imaging (circles) and 
2D images (squares) for a suspension with salt concentration 
[KCl] = 0.01 M. Note that the experimental second layer cor¬ 
responds to the sum of the second and third layers in the 
MC simulations. The layer boundaries in the numerical cal¬ 
culations are set at zi = 0.9do, and Z 2 = 1.8do and in the 
experiments are obtained from Gaussian fitting (see Fig. 


pled quasi-2D layers in thermodynamic equilibrium with 
respect to particle exchange. 

The equilibrium condition for layers i and z -I- 1 is 

Pi -I- mgzi = pi+i + mgzi+i, (24) 

where pi is the chemical potential of layer z, and Zi is 
its position. In our model, pi is approximated as the 
chemical potential of a 2D hard-disk fluid of area fraction 
(j)i. All disk diameters are equal to the sphere diameter 
d, which corresponds to a layer of spheres with the same 
vertical position z. 

In the low area-fraction limit, the chemical potential 
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Figure 11. Low-density limit of the near-wall particle dis¬ 
tribution for a monodisperse system (solid line) and polydis- 
perse systems with standard deviation of particle diameter 
(j/do = 0.1 (dashed), 0.15 (dotted) 0.20 (dash-dotted), and 
0.25 (long-dashed). 


of a hard-disk fluid is 


fi^ = kBT\ncl)^ + C{T), (25) 

where C{T) depends only on the temperature T. Ac¬ 
cording to the equilibrium condition (1241) and equation 
of state (1251) . we thus have 


(jii+i = 1 = 1 , 2 ,... (26) 

with the ratio r given by the Boltzmann factor 

r = e-'^/', (27) 

where I is defined by Eq. o and A = Zi+i — Zi. We 
assume that the layer separation A is independent of i. 
For finite area fractions, relation (|M1) is replaced with 

4)i+i = r{4>i)4>i, i = (28) 


where the layer occupation ratio r depends on the area 
fraction in the adjacent layers. The factor r{(j)) is deter¬ 
mined from the equilibrium condition (12411 with the help 
of the Gibbs-Duhem relation 

Trr/^ 

d^ = —(29) 

where p is the lateral 2D pressure within the layer. Com¬ 
bining (IMl) . (1^ . and (1^ yields 


where 


dr p'(,4>) 

d4> \_p'(r(j)) 



(30) 


(31) 



Figure 12. Occupation fraction fi for the first (black), second 
(blue) and third (green) layer vs the total area fraction (j) for a 
monodisperse HS suspension. Theoretical results (solid lines); 
MC simulations (symbols). 


The differential equation (1501) is solved for r = r{(j)) with 
the boundary condition (l27l) at (j) = 0. Occupation frac¬ 
tions fi = 4>i/(j) are then determined by iteration, apply¬ 
ing Eq. (I05|l and the relation 


OO 

= (32) 

i=l 

We have solved Eq. (l30ll and determined the occupa¬ 
tion fractions fi using the scaled-particle-theory equation 
of state for hard disks [5l|, 


TTcPp 4> 

4 ^ = 


(33) 


The results of our calculations are presented in Fig. [12] 
for a HS system with the same value of the sedimenta¬ 
tion length (HU as in our MC simulations. Based on 
the separation between the first and second peak of the 
suspension density profile shown in Fig. |5Ka), the calcu¬ 
lations were performed for A/d = 1. 

The theoretical results in Fig. [T5| are compared with 
the MC simulations of a monodisperse HS suspension 
with the boundaries between the layers set to zi = d and 
Z 2 = 2d, consistent with the peak positions. The agree¬ 
ment between our simple theory and simulations is quite 
good. A similar agreement was obtained for other values 
of the dimensional parameter I/d (results not shown). 

The layer boundaries used in Sec. IIVBI and IV Bl to 
compare the MC results with experiments differ from the 
boundaries used in the above model by approximately 
10%. Due to the observed deviation between the mea¬ 
sured and simulated particle distributions (see Fig. jH]), it 
is not possible to define the layer boundaries in a unique, 
equivalent way for the experimental and simulated sys¬ 
tems. Therefore, the layer boundaries Z\ = 0.9d and 
Z 2 = l-8d used in Sec. IIV Bl and IV Bi were chosen based on 
the comparison between the experimental and numerical 
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Figure 13. Short-time self-diffusion coefficient Ds, normalized 
by the Stokes-Einstein diffusion coefficient Do, as a function 
of the total area fraction 0; circles (squares) correspond to the 
first (second) layer, and triangles to the effective Ds calculated 
by the weighted average of the self-diffusion coefficients for the 
two individual layers, (a) Suspension with no added salt, (b) 
Suspension with salt concentration [KCl] = 0.01 M. 

results for the occupation fractions and self-diffusivities 
in particle layers. 


V. PARTICLE DYNAMICS 
A. Experimental results 

The short time self-diffusion coefficient in the x-y 
plane, is determined for different total area frac¬ 
tions of the sedimented particles by extracting the mean 
square displacement (l22ll from 2D epifluorescent images 
of the first and second particle layer. The mean-square 
displacement is measured over a time interval r that is 
small compared to the structural relaxation time of the 
suspension, to ensure that the measurements yield the 
short-time self-diffusion coefficient. 

The results are shown in Fig. [13] for suspensions with 
salt concentration [KCl] = 0.01 M and salt-free suspen¬ 
sions with [KCl] = OM. The self-diffusion coefficient is 
expected to decrease as the particle concentration in¬ 
creases; indeed, we observe this decrease for both salt 
concentrations and in both layers, for (f) < 0.4. In the 
case of [KCl] = 0.01 M, corresponding to A = 5 nm, the 
particles can get much closer to the cell floor, which in 


turn results in lower values of the self-diffusion coefficient 
compared to suspensions with [KCl] = 0 M. 

Using a linear fit to the values of Ds/Dq for the 
low area fractions, where there is no observable second 
layer, we can extrapolate to (p = 0 and extract the self- 
diffusivity of a single particle. The extrapolated results 
agree well with the measurements at very low concentra¬ 
tions (j) < 0.003, as discussed in Sec. IIV A II 

From the known occupation fractions fi and /2 for 
each (j) we can weigh the contribution of each layer to 
the total self-diffusivity, and construct an effective Dg of 
the whole suspension (Fig. [12]). As expected, for p < 0.4 
the effective self-diffusion coefficient Dg decreases as cp is 
increased in both salt concentrations. For larger p we ob¬ 
serve a flattening of Dg, which clearly indicates that the 
second layer becomes dominant in those area fractions. 
This observation is supported also by the saturation of 
at ^ > 0.45 [Fig. [3(b)] . 

B. Numerical simulations 

The results of our numerical simulations for the short- 
time lateral self-diffusion coefficient Dg in a HS system 
are presented in Fig. [TT] for a monodisperse suspension 
and for polydisperse suspensions with a/do = 0.1 and 
0.15. Figure [TS] shows the corresponding results for a 
system of monodisperse hard spheres with particle-wall 
and particle-particle electrostatic repulsion ([S]) and 
The results depicted in Fig. [TJ] indicate that for moder¬ 
ately polydisperse suspensions (in the range correspond- 



Figure 14. Normalized short-time self-diffusion coefficient 
Ds/Dq, as a function of the area fraction p for a suspen¬ 
sion with salt concentration [KCl] = 0.01 M. The main panel 
shows Ds averaged over the whole system, and the inset shows 
Ds for the first (bottom, red) and second (top, blue) particle 
layer. Experimental results (solid circles); simulation results 
(open symbols) for a monodisperse system (circles) and poly¬ 
disperse systems with the standard deviation of the particle 
diameter a = O.ldo (triangles) and 0.15do (squares). Note 
that at low area fractions the triangles overlap with the solid 
circles. The lines are a guide for the eye. 
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0 0.2 0.4 0.6 

Figure 15. Normalized short-time self-diffusion coefficient 
Ds/Dq, as a function of the area fraction ifi for a suspension 
with no salt. Symbols are the same as in Fig. 1141 Results are 
shown only for monodisperse suspensions. 



Figure 16. Normalized short-time self-diffusion coefficient 
Ds/Dq in the low-area-fraction limit 0 = 0 (for the system 
with salt) as a function of the suspension polydispersity. 


iiig to the polydispersity of silica particles used in the 
experiments), the self-diffusion coefficient is only moder¬ 
ately dependent on cr/do. For larger values of the vari¬ 
ance of particle diameters, the normalized self-diffusion 
coefficient Ds/Dq significantly increases with the degree 
of the polydispersity, because the mobility is dominated 
by small particles. This increase is illustrated in Fig. [TOl 
for a suspension in the low-area-fraction limit 0 = 0. 

The results of our hydrodynamic calculations for a 
HS suspension and for a suspension with screened elec¬ 
trostatic repulsion are compared with experimental re¬ 
sults for suspensions with [KCl] = 0.01 M (Fig. [T3]) and 
[KCl] = OM (Fig. [15]). For the system with salt, the 
measured values are slightly closer to their numerical 
analogs than in the absence of salt. Summarizing, the 
experimental and numerical results agree well both for 
the overall self-diffusivity and for the self-diffusivity in 
individual particle layers. 


VI. DISCUSSION 

In this paper we have studied in detail the structure 
and dynamics of quasi-2D colloidal suspensions near a 
wall, comparing experiment and theory. Our central re¬ 
sult is a rather sharp formation of a distinct second layer 
at an area fraction of 0 ^ 0.3. This value is much lower 
than the area fraction required for close-packing or other 
2D structural changes such as the formation of hexatic 
or crystalline order. One important consequence of this 
result concerns the apparent self-diffusion of the particles 
in the suspension and its dependence on particle density. 
Due to the higher mobility of the particles in the ele¬ 
vated layer, the effective diffusivity is higher and levels 
off as particle density increases. The experimentally ob¬ 
served behavior could be interpreted incorrectly if one is 
unaware of the layering (or stratifying) effect. 

We find good agreement between experimental and 
simulation results for the occupation fractions of the 
first and second layers and for the lateral self-diffusivity 
(both for the entire suspension and in the individual lay¬ 
ers). However, we also find an unexpected discrepancy 
in the position and the height of the second peak in the 
near-wall particle distribution. While the source of this 
discrepancy is unknown, one possibility, related to op¬ 
tical aberrations, is suggested in Appendix [^ On the 
other hand, the difference between theory and experi¬ 
ment might also be a result of an actual physical effect, 
such as more complicated electrostatic interactions set¬ 
ting in at higher layer densities. 

Another new insight put forth in this study is the sig¬ 
nificant effect that polydispersity has on the occupation 
and composition of layers close to the bottom wall, even 
in the case of a relatively small dispersion of particle sizes. 
The effect of polydispersity is evident already at low den¬ 
sities, since the smaller and larger particles segregate into 
the upper and lower layers, respectively. We expect the 
phenomena described here to be quite general and to be 
manifested in any such system where the sedimentation 
length I is of the order of the particle diameter. This con¬ 
clusion is supported by the appearance of the phenomena 
both in experiments and in Monte-Carlo simulations. 

An important outcome of this paper is the construction 
of a very simple theoretical quasi-2D model of the layered 
microstructure in thermodynamic equilibrium. Such sys¬ 
tems have been analyzed earlier using density-functional 
theory [^ . but our theoretical model is much simpler 
and easier to apply. We have demonstrated that the 
model approximates well the experimental and numer¬ 
ical results for the system studied in this work. 

We conclude with three open issues. Layering phenom¬ 
ena near a wall are well documented in 3D suspensions as 
well [in. An interesting question is whether this per¬ 
turbation to the 3D pair correlation function could be 
fundamentally related to the sequential layering reported 
here. The structural features near the wall should also 
affect two- and many-particle dynamics in the quasi-2D 
suspensions, which can be characterized by two-point mi- 
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crorheology. Finally, taking a more detailed account of 
interparticle forces such as strong electrostatic interac¬ 
tions may hopefully provide deeper understanding of the 
effects observed in this work. 
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Appendix A 

We present a simple model to support a hypothesis 
that the discrepancy between the measured and calcu¬ 
lated near-wall particle distributions stems from optical 
aberration caused by nonuniform optical properties of the 
suspension in the near-wall region. We assume that such 
aberration produces a nonlinear rescaling of the coordi¬ 
nate 


z = z{z), (Al) 

where z is the actual and z is the measured particle po¬ 
sition. The rescaling (EH) results in the corresponding 
transformation of the particle density 

p{i)=p{z)^. (A2) 

a z 

To demonstrate that a distortion (ED can produce the 
observed shift and change of height of the features of the 
distribution p, we consider an ad hoc distortion model 
with the transformation between the measured and ac¬ 
tual vertical coordinates given by the equations 


dz 

dz 


and 


1 ; 

1 - (1 - 5) 
b] 


z — Zl 

Z2 - Zl’ 


Z < Zl, 

Zl < z < Z 2 , (A3a) 

Z2 < -Z, 


b = a(j), (A3b) 

where zi, Z 2 , and a are the model parameters. The trans¬ 
formation (IA3I) describes position-dependent coordinate 
contraction with the amplitude gradually increasing in 
the region zi < z < Z 2 (the region where the second peak 



z/dg 



Figure 17. The effect of the model distortion transformation 
(IA3I) on the particle distribution p{z) in a HS suspension with 
the standard deviation of the particle distribution a/do = 
0.2. (a) A comparison of the MC result for p(z) (solid line) 
with the transformed distribution (lA^ (dashed line) at the 
area fraction = 0.54; (b) the transformed distribution for 
different area fractions (as labeled). The parameters of the 
transformation (IA3l) are zi/d = 0.55, Z 2 /d — 1.2, and a = 2.5. 


occurs according to the experimental data). The overall 
deviation of the Jacobian (IA3I) from unity is proportional 
to the area fraction of the suspension layer. 

Figure [TTl al compares the distorted distribution (jA2l) 
with the corresponding untransformed distribution p{z) 
obtained from MC simulations of a HS suspension at the 
area fraction (f) = 0.54. Figure [TTIb) presents the dis¬ 
torted distribution for the set of area fractions for which 
experimental results are depicted in Fig. [SDa). The pa¬ 
rameter values of the transformation (IA3I) are given in 
the hgure caption. 

The results show that the coordinate transformation 
(IA3al) shifts the position of the second particle layer to 
the left and produces a corresponding enhancement of 
the peak of particle distribution, similar to the experi¬ 
mentally observed features of the distributions depicted 
in Figs.EKa) andEDb). Thus our calculations provide in¬ 
direct support to our optical-distortion hypothesis. The 
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distortion hypothesis can also explain why the measured 
and calculated occupation fractions and self-diffusivities 
of the particles in the top layer agree well (see Figs. flOl 
and [13]), in spite of the fact that the observed and calcu¬ 
lated positions of the layer differ significantly. 

It is an open question what the source of the distor¬ 
tion (EH) might be. Since the suspension is imaged from 


above in our confocal-inicroscopy system, we hypothesize 
that reflection of laser light from the first (bottom) parti¬ 
cle layer results in stray illumination of the second layer, 
producing distorted particle height measurements. The 
optical distortion hypothesis can be verified by exper¬ 
iments using refractive-index matched suspensions, but 
such investigations are beyond the scope of the present 
study. 
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